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Abstract 

Stability arguments suggest that the Kaluza-Klein (KK) internal scalar 
field, should be coupled to some external fields. An external bulk real scalar 
field, ip, minimally coupled to gravity is proved to be satisfactory. At low 
temperature, the coupling of ip to the electromagnetic (EM) field allows <I> to 
■^j- ■ be much stronger coupled to the EM field than in the genuine five dimensional 

KK theory. It is shown that the coupling of <I> to the geomagnetic field may 
explain the observed dispersion in laboratory measurements of the (effective) 
gravitational constant. The analysis takes into account the spatial variations 
of the geomagnetic field. Except the high PTB value, the predictions are 
^ | found in good agreement with all of the experimental data. 



PACS numbers: 04.50. +h, 95.30.Sf, 06.20.Jr, 04.80.Cc 

1 Introduction 

Although the methods and techniques have been greatly improved since the late 
nineteenth century, the precision on the measurement of the gravitational constant, 
G, is still the less accurate in comparison with the other fundamental constants 
of nature 0. Moreover, given the relative uncertainties of most of the individual 
experiments (reaching about 10 -4 for the most precise measurements), they show 
an incompatibility which leads to an overall precision of only about 1 part in 10 3 [[|. 
Thus the current status of the G terrestrial measurements (see HJ) implies either an 
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unknown source of errors (not taken into account in the published uncertainties), 
or some new physics [|]]. In the latter spirit, many theories have been proposed as 
candidates for the unification of physics. As such, they involve a coupling between 
gravitation and electromagnetism (hereafter GE coupling), as well as with other 
fields present. 

The Kaluza-Klein theories [||, || have been among the first attempts to unify 
electromagnetism and gravitation. Although they are disqualified in their original 
form, they constitute the prototype for many present theories (in particular, those 
involving extra-dimensions). In the laboratory conditions for G-measurements, fun- 
damental theories as well as KK theories are well described as an effective theory, 
under the form of the Einstein and Maxwell equations conveniently modified, as 
described below. 

This paper explores the possibility that the discrepancy between the results of 
the G-measurements is the effect of the GE coupling, by comparing the available 
data to the predictions of such an effective theory. We do not pretend that the 
effective theory adopted here (KK0, see below) is the only possibility. Rather, we 
consider that it provides the simplest opportunity to confront the idea of a GE 
coupling with real data. The success of the fit suggests that this may be the real 
explanation. 

The simplest theories accounting for a GE coupling are those of Kaluza-Klein 
H U, or slight modifications of them. Their use corresponds to the most economical 
way (the minimum of hypotheses) to test the hypothesis of GE coupling. As it is well 
known, such theories are effectively well described by Einstein equations (with their 
correct Newtonian limit), where the newtonian gravitational constant G is replaced 
by G e ff, given below: this effective gravitational constant depends of the fifteenth 
degree of freedom, (744, of the (5-dimensional) bulk metric, which plays the role of a 
(four dimensional) scalar field $ (#44 = — $ 2 in the Jordan-Fierz frame). 
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The genuine five dimensional Kaluza-Klein theories being subject to instabilities 
(IS H> H)i various authors ([|H], |ll|, [12|]) have suggested a more acceptable version 
which includes an additional stabilizing external bulk field: here we adopt the min- 
imal hypothesis of a scalar field ip minimally coupled to gravity. In this theory 
(hereafter KKip), G e ff varies with the electromgnetic field, and thus in spacetime. 
However, only variations with respect to the cosmic distance or time have been in- 
vestigated in the literature hitherto. But, since the geomagnetic field varies, with 
latitude and longitude, and thus at the different sites of the G measurements, The 
GE coupling implies that the experiments in fact measure the distinct corresponding 
values of G e //, rather than an unique value of G. 

Note also that this theory predicts a variation of the effective fine structure 
constant a with the gravitational field, and thus with the cosmological time. In 
a companion paper (fl3]|), we compare (with success) the predicted evolution with 
astrophysical data concerning the distant quasars ( [JT^j, |Hi[|). 

In section 2 we recall the definition of the effective coupling constants, and the 
effective Maxwell-Einstein equations, issued from the five dimensional compactified 
KK theory stabilized by a minimally coupled bulk scalar field. In section 3, we 
calculate the vacuum solutions on Earth, in the weak field limit, taking into account 
the geomagnetic field. In section 4, we confront the predicted values of G e ff with 
the laboratory measurements. In section 5, we discuss the consistency of our results 
with respect to the orbital motion of the LAGEOS satellite, the Moon and planets 
of the solar system, as well as the binary pulsar PSR1913 + 16. 

2 Theoretical background 



An argument initially from Landau and Lifshitz |16j may be applied to the pure 
Kaluza-Klein (KK) action (|12[) : the negative sign of the kinetic term of the KK 
internal scalar field leads to inescapable instability. Stabilization may however be 



obtained if an external field is present (J7|, ||, |T2[), and we assume here a version 
KKip of the KK theory which includes an external bulk scalar field minimally coupled 
to gravity. After dimensional reduction (a = 0,1,2,3), this bulk field reduces to 
a four dimensional scalar field ip = ip(x a ) and, in the Jordan-Fierz frame, the low 
energy effective action takes the form (up to a total divergence) 



c 4 $ 1 a c 4 d $<9 Q $ 



+ / v /T 9*[ia«'/' - f - J</>]<i 4 i, (1) 

where v4 Q is the potential 4- vector of the electromagnetic field, F a p = d a Ap — dp A a 
the electromagnetic field strength tensor, U the self-interaction potential of ip and 



J its source term. Following Lichnerowicz 1/], we interpret the quantity 



Geff = f (2) 

of the Einstein-Hilbert term, and the factor 

e 0eff =e $ 3 (3) 

of the Maxwell term respectively as the effective gravitational "constant" and the 
effective vacuum dielectric permittivity. The effective vacuum magnetic permeability 
reads /io e // = Ho/& 3 , so that the velocity of light in vacuum remains a true universal 
constant. Both terms depend solely on the local (for local physics) or global (at 
cosmological scale) value of the KK scalar field $, assumed to be positive defined. 

The source term of the ip-G.eld, J, includes the contributions of the ordinary 
matter (other than the scalar fields ip and $), of the electromagnetic field and of the 
internal scalar field For each, the coupling is defined by a function (temperature 
dependent, as for the potential U) fx = fxijp, $), where the subscript X stands for 
"matter", "EM" and In order to recover the Einstein-Maxwell equations in the 
weak fields limit, these three functions are subject to the conditions: /bm^ , 1) = 



fmatter{v , 1) = /$ (i> , 1) = 0, where v denotes the vacuum expectation value (VEV) 
of the -0-field. 

The contributions of matter and $ are proportional to the traces of their re- 
spective energy-momentum tensors. Since the energy-momentum tensor of the elec- 
tromagnetic field is traceless, a contribution of the form e Jem F a p F al3 accounts 
for the coupling with it. The fit of our model to the data (see below) shows that 
9 g S ^ 1 (v , l)v> AttG/c 4 , as it can be expected near the vacuum at low temperature. 
Thus, we will not take the latter term into account. However, we may suspect that 
dg,$ A/ (v , 1) v < ArcG/c 4 at high temperature, which may have consequences in some 
astrophysical conditions (see below). 

Applying the least action principle to the action ([[]) yields : 



the generalized Einstein equations 

p 1 F> n - 8nGe ff [t( £M ) _i_ t<*' _i_ tWi ' n 

tt a p--Kg a p- — — A — [i af3 + i aj3 + l a/3 J, (J) 



where 



T ip M) - £ 0eff i.~ F 2 F M + 7 F -yS F1& 9a/3 ), (5) 



^ 4 



S 8^G (VaV/3$ -9apW$) ( 6 ) 



and 
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the generalized Maxwell equations 

V 7 F a p + V/3 F 7Q + V a Fpy = 

and 



V" ( £ e// Fa/3 ) = 0, (9) 
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and the scalar fields equations 



and 

„ ™^ 4vrG 
V„ V $ = - 



(9 7 1 

(11) 

where the symbol V„ stands for the Riemannian covariant derivative. Clearly, 
T^ M ^ and T^g define respectively an effective energy-momentum tensor for 
the electromagnetic field in the presence of the KK scalar field and an effective 
energy-momentum tensor for the latter itself. Relations are formally the 
same as the Einstein-Maxwell ones, but with the additional contribution of 
the KK scalar as a matter source and the replacement of G and Eq by their 
respective effective values. 

3 Vacuum solutions in the presence of a dipolar 
magnetic field 

Since we are in weak field conditions (we look for small deviations from Newtonian 
physics), we only keep first order terms. Thus, we neglect the excitations of $ and 
ip with respect to their respective VEV's 1 and v. Also, the energy density of the 
■0-field must be lower than that of the magnetic field. 

Let us study the spatial variation of xp out of the fields' source, but in presence of 

a static dipolar magnetic field, B = B(r,(p,9). We denote r, (p and 9 respectively 

the radius from the centre, the azimuth angle and the colatitude. Thus, writing 

ip = ip(r, if, 9) and $ = $(r, (p, 9), and taking into account that ) = (definition 

of the VEV), equations (0) and (JTTJ) simplify respectively as 

l,d dip 1 d dip 1 0V, 9f EM B 2 

r 1 Or or sin 09 09 sin 9 dip 2 dtp fj, 



and 

1 d , 2 d® 1 d d$ 1 <9 2 $ df EM B 2 

~9 "T" ) H : — J. ( Sin 9 -p~r ) H T~2j. = — 2 — {v , l)v — , (13) 

v dv dr sin 9 39 89 sin 2 <9y? 2 9$ /x 

where we have dropped the pure gravitational constant 47rG/c 4 with respect to 
v 9 Jem I d^, as indicated above. 

We consider a dipolar magnetic field: B — V V. For our purpose, it is sufficient 
to limit the expansion of the scalar potential, V, to the terms of the Legendre func- 
tion of degree one (n = 1) and order one (m = 1). Hence V = (a 3 /r 2 ) [g® cos9 + 
g\ sm.9 cos ip + h\ si~n.9 siny], where g®, g\ and h\ are the relevant Gauss coeffi- 



cients, a is the Earth's radius and M = y(9i) 2 + (fii) 2 + (M) 2 denotes its 



magnetic moment (see |18], jI5f). Setting cosy?i = — g\j ' \i \g\) 2 + (h\) 2 , simp 1 



h\/\J (g\) 2 + (h\) 2 and tan A = g\j ' \J (g\) 2 + (h\) 2 , the solutions of equation (|13|) 
then reads 

$ = 1 - k{r)x(9,cp), (14) 



where we have set 



and 



fc ( r ) = — ^-(^i)ni — n) (is) 



yU <9$ 47r r 



a;(6»,v9) = cos 2 6> + e(9,<p) 
= cos 2 6* + cot 2 A sin 2 cos 2 ( <p + ipi) — cot A sin 29 cos ( <p + ipi ) (16) 
and similarly for ^ by making the substitution — > — in relation (|T5|). 



Thence, one derives the expression of G e ff(r,9,(p) by inserting the solution ( |i4|) 
above in relation (|2|). Thus, since (v , 1) t> > and the variable x turns out to 
be positive at any position in space, it follows that the effective gravitational constant 
G e ff will always be greater than the true gravitational constant, G. Whence the 
prediction of an upward bias in the laboratory measurements of G. 
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4 Comparison with laboratory measurements 

Because of various uncontrolled systematic errors, the data published by the different 
laboratories have different precisions. In the following, we test two hypotheses with 
respect to these results: HO = Hypothesis of a constant G (i> — n — 1) and HI = 
Hypothesis of an effective G {y = n — 2). Here v denotes the number of degrees of 
freedom, n is the number of data points, and we have 1 or 2 parameters in the fit. 

There are presently almost 45 results of measurements G published since 1942 
(see e.g., |§, Table 2, pp. 168 and 169). We exclude from the present study the mine 
measurements because of the too numerous uncontrolled systematic biases involved. 
The " accepted " values are presently G = 6.67259 ± 0.00085 lO" 11 (CODATA 86, 
H) andG = 6.670±0.010 10~ n (CODATA 2000, |2lJ) in MKS unit. A fitting of the 
45 data with these values give respectively xl — 145.17 and xl — 213.25 (xl — X 2 / v -> 
where v denotes the degrees of freedom). If we forget the accepted value and try a 
best fit, assuming an arbitrary constant value of G, we obtain G = 6.6741 10 -11 SI 
with xl = 127.04. 

The more discordant laboratory measurement (high PTB value |ITJ) is contro- 
versial. If we discard it, the previous fits lead to xl = 11128 (CODATA 86), 
Xl = 62.498 (CODATA 2000) and xl = 2.341 (free value). Since it seems now 
certain that this high PTB value |4l| suffers from some systematic error (see p6[ , 



for more details), we remove it for our analysis (if we keep it, our model is still 
more favored). Note that the strongest contributions to the xl then come from the 
BIPM 2001 [|36| and the HUST [F25J measurements. Thus, unless there are some 
experimental systematic errors presently not understood, these experiments do not 
measure the same quantity. In the framework of the theory proposed here, they 
measure G e //- Because of the GE coupling, G e // should depend on the geomagnetic 
field at the laboratory position. 
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First, we select a subset of results with a good precision and which does not 
suffer any significant systematic error: we include only those points with relative 
uncertainty SGi a b/Gi a b < 1CT 3 . Also, we demand a short measuring time (At < 
200 s), to minimize the possible biases due to the time variations of the geomagnetic 
field. This gives the sample SI, with 17 points, shown in figure 1. We fit these 
data using the IGRF 2000 Gauss coefficients^ (<# = -0.31543, g\ = -0.02298 and 
h\ = 0.05922 in Gauss). We obtain (Figure 1, Table 2) 

— = (0.149933 ±0.000017) - ( 0.0001514 ± 0.0000262 ) x{L, I), (17) 

10 G e ff 

in MKS units, with xl — 1-327. This gives 

G = ( 6.6696 ± 0.0008 ) 10~ n m 3 kg" 1 s~ 2 , (18) 

that we retain as a true gravitational constant. The relative uncertainty is only 1 
part in 10 4 : the major part of the differences between the laboratory measurements 
was generated by the predicted variation of G e ff with the magnetic field. Further, 
adjusting to the same set the mean value M = 7.87 10 22 A m 2 (in the time interval 
spanning from 1942 to 2001), it follows 

(v , 1) v = ( 5.44 ± 0.66 ) 10~ 6 fm Tev~\ (19) 

that we retain too. We observe also that the HUST value (the lowest most precise 
measured value of G) is perfectly fitted: it differs from other values because of the 
proximity of this laboratory to the equator. 

Then, we fit the whole sample (excluding the PTB value, as stated above). We 

lr The IGRF coefficients are given for time intervals of five years. Hence, for a more precise 
fitting, one should use the most suitable IGRF coefficients for a given laboratory value according 
to the years the measurements were carried out. Of course, the necessity of doing so depends on 
the precision reached in the laboratory G measurements. Nevertheless, we have also computed 
the variable x(L, I) using the IGRF 1965 Gauss coefficients (g\ = - 0.30339, g\ = - 0.02123 and 
h\ = 0.05758 in Gauss). However, no significant changes to our conclusions were found. 
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obtain the set of the 44 measurements given in figure 2: the best fit leads to 



— — — — = ( 0.149929 ± 0.000017 ) - ( 0.0001509 ± 0.0000252 ) x(L, I) (20) 
10 G eff 

in MKS units. It gives xl — 1-669, to be compared to xl — 2.255 for the best fit 
assuming a constant G. These results are summarized in Table 2. In order to check 
the relevance of our result (which involves two free parameters rather than one), we 
apply the F test (Fisher law). This yields F x = = 16.09, which indicates that, 
independently of the number of parameters, our fit is better with a significance level 



gre ater than 99.9% 



Location [reference] 


Latitude (°) 


Longitude (°) 


G lab (10- 11 m 3 kg- 1 s- 2 ) 


Lower Hutt (MSL) f§, || 


-41.2 


174.9 


6.6742 ± 0.0007 
6.6746 ± 0.0010 


Wuhan (HUST) || 


30.6 


106.88 


6.6699 ± 0.0007 


Los Alamos |p6(| 


35.88 


-106.38 


6.6740 ± 0.0007 


Gaithersburg (NBS) |27|, 




38.9 


-77.02 


6.6726 ± 0.0005 
6.6720 ± 0.0041 


Boulder (JILA) ]2g 


40 


-105.27 


6.6873 ± 0.0094 


Gigerwald lake |0|, 


46.917 


9.4 


6.669 ± 0.005 (at 112 m) 
6.678 ± 0.007 (at 88 m) 
6.6700 ± 0.0054 


Zurich H, || 


47.4 


8.53 


6.6754 ± 0.0005 ± 0.0015 
6.6749 ± 0.0014 


Budapest B4 


47.5 


19.07 


6.670 ± 0.008 


Seattle || 


47.63 


- 122.33 


6.674215 ± 0.000092 


Sevres (BIPM) ||, 


37 




48.8 


2.13 


6.67559 ± 0.00027 
6.683 ± 0.011 


Fribourg |38 


46.8 


7.15 


6.6704 ± 0.0048 (Oct. 84) 
6.6735 ± 0.0068 (Nov. 84) 
6.6740 ± 0.0053 (Dec. 84) 
6.6722 ± 0.0051 (Feb. 85) 


Magny-les-Hameaux 


3S] 


49 


2 


6.673 ± 0.003 


Wuppertal || 


51.27 


7.15 


6.6735 ± 0.0011 ± 0.0026 


Braunschweig (PTB) gT], || 


52.28 


10.53 


6.71540 ± 0.00056 
6.667 ± 0.005 


Moscow 43, 44 




55.1 


38.85 


6.6729 ± 0.0005 
6.6745 ± 0.0008 


Dye 3, Greenland [|45J| 


65.19 


-43.82 


6.6726 ± 0.0027 


Lake Brasimone |46| 


43.75 


11.58 


6.688 ± 0.011 
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Table 1 : Results of the most precise laboratory measurements of G published 
during the last sixty years and location of the laboratories. 



Sample 


HO 


HI 


SI 

17 points 
(Fig.l) 


Xl = 3.607 (best fit) 
Xl = 21.523 (mean of CODATA 86) 
Xl = 141.46 (mean of CODATA 2000) 


Xl = 1.327 


Whole 
[23] - 47] 
44 points 

(Fig.2) 


Xl = 2.255 (best fit) 
Xl = 11.128 (mean of CODATA 86) 
Xl = 62.498 (mean of CODATA 2000) 


Xl = 1.669 



Table 2 : Reduced x 2 f° r the two different hypothesis HO (Hypothesis of a con- 
stant G) and HI ((Hypothesis of an effective G), and different samples SI and 
whole (except the high PTB value [O], see text). 



5 Discussion and conclusion 

We have not taken into account the temporal variation of the geomagnetic field. 
All the data considered in this paper are averaged on time. Nevertheless, for those 
which would not be averaged on time, one may expect small time variations of 
G e ff both with the Sq and with the L field disturbances of the geomagnetic field. 
Now, periodic variations of the gravitational "constant" with the lunar or diurnal 
period have yet been pointed out in the literature (see j|8], |26|, (43|). Although 
it is presently believed that they are related to tides, the explanation could be this 
temporal variation. We notice that the G measurements of |38j are consistent with 
an annual variation. 



Besides, a recent study j|9| shows that helio seismology seems to favor a low 
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(3 




Figure 1: Laboratory measurements with relative uncertainty s ^ lab < 10 3 and 



measuring time At < 200 s (sample SI, 17 points |3|, H - IE 0, @, H 



- HO], [^3] - [0). The line indicates the best fit G\ a \, versus the mixed variable 
(xl = 1-327). Assuming a constant G would yield a bad fit to the data (xl — 3.607) 
mostly because of the HUST value. 



x 



6,7400 
^ 6,7200 

CO 

i? 6,7000 

CD 

£ 

t 6,6800 
o 

O" 6,6600 
6,6400 



PTB 95 



o 




0,2 0,4 0,6 

x(L,l) 



0,8 



Figure 2: versus x (whole sample plus the PTB 95 value, 45 points 
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value of G, close to the HUST value. This is in accordance with our predictions (see 
section 3 and relations fll7D and (f[8| )), since one expects the effective coupling of 
the $-field to the EM field to decrease towards its lowest value as the temperature 
increases. Hence, one gets a good agreement between our prediction ([Tj^) and the 
helioseismic data, on account of the high temperature met in the core of the Sun. 

To test our proposal of a spatial dependence of G e ff, we call for a complete cov- 
ering of the Earth (in particular in the south hemisphere) and at different latitudes. 
On the basis the current available data, the present study allows to make some pre- 
dictions. First, according to the authors fl4*9 |, the Sudbury Neutrino Observatory 
(SNO: L = 46° 29' N and I = 80° 59' O) Q should constitute a promising means 
of determining G with a quite good accuracy. Hence, since the measured quantity 
is actually the effective gravitational constant, the value predicted at SNO by the 
relation flTTD is G S no = ( 6.6742 ± 0.0009 ) 10~ n m 3 kg" 1 s~ 2 , that is practically 
the same as at the MSL P5[ . In addition, it would be of interest to test also our 
predicted value at the surface of the poles, that is: G po i es = ( 6.6763 ± 0.0010 ) 10~ n 
in MKS units. At the equator, one should get a dependence with respect to the 
longitude (since the magnetic and geographic poles do not coincide). 

At the orbit of the Moon, the relative deviation, (G e // — G)/G, to the true 
gravitational constant is predicted to be as small as 1.7 10~ 13 which is consistent 
with lunar laser ranging (see, f5"I|| ). In the Solar System, the relevant magnetic 
field is the dipolar field of the quiet Sun. Since the coupling constant of the $- 
field to the EM field is much weaker within the Sun than on Earth, one finds at 
the orbital radius of Mercury [G e jf — G)/G < 10~ 6 and much smaller beyond, 
decreasing as 1/r 4 . Thus, at the orbital radius of Neptune (G e // — G)/G drops to 
less than 10 -12 . Taking into account the overall planetary constraints on GM Q \f>2\ . 
the accord between the proposed model and observation is still acceptable. 

For artificial satellites in quasi-circular orbits, the appropriate quantity for com- 
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parison with observational data is (G e ff — G')/G' when the analysis is based only on 
the satellite data (see appendices A and B), where we have set G' = G [ 1 + ^p- ( 13 + 
| ( 2 ) 2 J 2 ) ] and J2 is the true Earth quadrupole moment coefficient. Referring to 
the orbital motion of the LAGEOS satellite, with semimajor axis equal to 12 270 km, 
inclination i = 109.94° and eccentricity e = 0.004, one finds the maximum devia- 
tions (G eff - G')/G' ~ 1.9 10- 8 (whereas (G eff - G)/G ~ (G' - G)/G ~ 2 lO" 5 ), 
consistent with the constraint | a \ < 10~ 5 — 5 10~ 8 on the Yukawa coupling constant 
|53| , figure 3.2a, p. 99 and section 6.7). 



a see 



We also looked for a possible relative deviation to the true gravitational constant 
induced by the strong magnetic fields of pulsars. We found this effect quite negli- 
gible, of the order 10~ 7 in the case of the binary pulsar PSR1913 + 16, for which 
observation yields G e // = Gn (l-OO^o^) The deviation is tiny because of the 
small radius of the pulsar and of its companion as compared to their respective 
orbital radius around the center of mass. 

Thus we conclude that, apart from systematic errors that need to be corrected 
(e.g., by applying some prescriptions like that pointed out by Kuroda |]55) for the 
swinging pendulum method), the only possibility to reconcile the published values 
of G is to consider a dependence on the latitude and longitude, of the type proposed 
here. In particular, if all present systematic errors could be removed in the future, 
we predict Gi ab = (6.6742 ± 0.0009) 10 _u m 3 kg" 1 at the PTB laboratory, that is 
the same value as that predicted at SNO and the current MSL. 

Up to now, a lot of attention had been paid to the dependence of G on cosmic 
time or radial distance only. But the dependence on latitude and longitude, that we 
examine here, has not been taken into account. More precise measurements (e.g., 
the SEE project |56[]) and further analyses taking into account higher harmonics 
and the various kind of changes of the geomagnetic field should bring more support 
to our claim. 



14 



6 Appendix A 

The generalized Einstein equations (f|) rewrite (including the contribution of the 
ordinary matter, T^) 

+ (T$ - \T^g aP )\. (21) 

Now in the weak fields and slow motion approximation, only the 00 components are 
relevant. Hence, the above equations reduce to 

A($ soo) = V$ . V In J=g + ^ [TT ] + ( Tff - \ T W goo ) + ( # - ~ T (m) *» ) ]• 

(22) 

The first term of the right hand side of the above equation is second order. Hence, 
neglecting the energy densities of the EM field and the V'-field with respect to the 
(ordinary) matter density, one gets 

A($Soo) = ^(T ( m) - \T^ g oo\ (23) 

Clearly, this yields the Newtonian potential divided by $, that is an effective po- 
tential where G^ is replaced by G e ff = G/<3>. 

7 Appendix B 

As a consequence, the equation of motion of a satellite should be sensitive to the 
effect of the effective G if present. Indeed, in the static gravitational field of a 
rotating body of mass M with angular velocity to and for r ~ 2a or more, one may 
merely replace in the first approximation G^ by G e ff in the Newtonian potential. 
Thus, d 2 r/dt 2 = — VV" = — (g r u r + geue), where u r = f/r, uq = du r /d9 and 

V = - ( ^-\\ - {-f hdcosH -\)\ - c 2 r(l - cos 2 9). (24) 
r z r 2 2 
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Hence, inserting relations ([L4| - |HJ) in equation ( pop above and expanding, one gets 
, r = ^^{l-3(^) 2 Ja^co^-i) + Mr)e(^,^)[l-3(^) 2 J 2 (^co^-i 

- - Mr) ( - f J 2 cos 4 6 } - oo 2 r ( 1 - cos 2 6 ) (25) 
2 r 

and analogously for gg, where 

jt = 3 \_L± 9 a 26 

1 + *fel[l 3 + |(fl)2 j 2 ] 

appears as the effective quadrupole moment coefficient of the central body at radius 
r. The quadratic cosine term provides an additional term to the effective zonal 
harmonic coefficient of order 4. Clearly, by interpreting the data from a single 
satellite in circular orbit, J' 2 will appears as a constant parameter. 
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